Backbending and 7 vibrations 
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We propose that the backbending phenomenon can be explained as a result of the disappearance 
of collective 7-vibrational mode in the rotating frame. Using a cranking+random phase approxi- 
mation approach for the modified Nilsson potential + monopole pairing forces, we show that this 
mechanism is responsible for the backbending in ^^^Dy and ^^**Er, and obtain a good agreement 
between theoretical and experimental results. 
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There is a general persuasion that the backbending is 
caused by the rotational aUgnment of angular momenta 
of a nucleon pair occupying a high-j intruder orbital near 
the Fermi surface. It is assumed that the alignment 
breaks a singlet Cooper pairing in this pair and leads 
to a sudden increase of the kinematical moment of in- 
ertia J^^^^ = I/Q along the yrast level sequence as a 
function of a rotational frequency (cf RefUl). Indeed, 
in many cases this single-particle (quasiparticle) mecha- 
nism is supported by microscopic analysis in terms of var- 
ious cranking Hartree-Fock-Bogoliubov calculations (cf 
0, 0, 3). It should be noted, however, that the role 
of vibrational (collective) excitations in the backbending 
has never been studied. 

It is well known that mean field description of finite 
Fermi systems could break spontaneously one of the sym- 
metries of the exact Hamiltonian, the so-called spon- 
taneously symmetry breaking (SSB) phenomenon. In 
fact, the mean field description of the backbending corre- 
sponds to the first-order phase transition. However, this 
concept is appropriate only in the limit of infinite number 
of particles. Obviously, for finite systems quantum fluc- 
tuations, beyond the mean field approach, are quite im- 
portant. The random phase approximation (RPA) being 
an efficient tool to study these quantum fluctuations (vi- 
brational and rotational excitations) provides also a con- 
sistent way to restore broken symmetries. Moreover, it 
separates the collective excitations associated with each 
broken symmetry as a spurious RPA mode and fixes 
the corresponding inertial parameter. In this paper we 
present the first self-consistent quantitative treatment of 
these ideas for rotating nuclei We demonstrate that 
the backbending in ^^^Dy and ^^^Er can be explained as 
a result of vanishing one of the quadrupole vibrational 
modes in the rotating frame at a critical rotational fre- 
quency. Consequently, collective motion associated with 
this mode should describe rotational states of the nonax- 
ial rotating system. 

The practical application of the RPA for nonsepara- 
ble effective forces such as the Gogny or Skyrme interac- 
tions in rotating nuclei requires too large configuration 
space and still is not-available. The RPA with separable 
multipole-multipole interaction based upon phenomeno- 



logical cranking Nilsson or Saxon- Woods potentials with 
pairing forces gives a sufficiently good description of low- 
lying collective excitations in rotating nuclei (cf |E0)- 
Following this approach, which hereafter is called CRPA, 
we start with the cranking shell model (CSM) Hamilto- 
nian in the form 



Hn = H- 



J2 KNr - 



(1) 



T = N,P 



The term H = Hj\j + i?add contains a Nilsson Hamilto- 
nian Hn with three different oscillator frequencies — 



MP, if 



[i = 1,2,3) that de- 



l-2/3y'^cos(7-f^) 
termine quadrupole deformation parameters 7 and /? (cf 
8J). The frequencies are subject to the volume conser- 
vation constraint LOxLOyi^z = i^oo (^"^oo = 41yl~3 MeV) 
that models the nuclear incompressibility. In the crank- 
ing model with the standard Nilsson potential the value 
of the moment of inertia is largely overestimated due to 
the presence of the velocity dependent / ^ term. This 
shortcoming can be overcome by introducing the addi- 
tional term i/add = I]f ^add(i) with 



hadd = rim cJooK 



(2) 



TOCJoo 2 



The term restores the local Galilean invariance of the 
Nilsson potential in the rotating frame and removes the 
spurious effects of the I ^ term (see for details Ref0) ■ 
Note that this basic recipe supersedes the fitting proce- 
dure of nuclear inertial properties used, for example, in 
Reftol 

The interaction is taken in a separable form 



^^int — — GrPr Pt 



1 



:K2 



(3) 



Here, r = neutrons or protons, P+ = c^c^ and c^, Cfc 
are creation and annihilation single-particle operators. 
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respectively. An index k is labelling a complete set of 
the oscillator quantum numbers (|fc) = \Nljm)) and the 
index k denotes the time-conjugated state We recall 
that the K quantum number (a projection of the angu- 
lar momentum on the quantization axis) is not conserved 
in rotating non-axially deformed systems. However, the 
CSM Hamiltonian adheres to the D2 spatial symmetry 
with respect to rotation by the angle tt around the ro- 
tational axis X. Consequently, all rotational states can 
be classified by the quantum number called signature 
a — exp(— iTra) leading to selection rules for the total 
angular momentum / = a+2n, n = 0, ±1, ±2 .... In par- 
ticular, in even-even nuclei the lowest rotational (yrast) 
band characterized by the positive signature quantum 
number cr = 1 (a = 0) consists of even spins only. 

The quadrupole operators Qm (™ = 0,1,2) are de- 
fined by 



^2+m+(CT+3)/2 
V/2(1 + 6rao) 



(4) 



The monopole interaction is defined by the positive sig- 
nature operator M^^^ = t^Yq. Single-particle matrix el- 
ements of any one-body Hermitian operator F {P, Q, M) 
are determined by the signature, time-reversal and Her- 
mitian conjugation properties of the operator (cf Ref|^. 
All multipoles are expressed in terms of the double- 
stretched coordinates qi = ^li, {<li = x,y,z). The ef- 
fective interaction restores the rotational invariance of 
the Hamiltonian H in the limit of the harmonic oscil- 
lator potential This is especially important in or- 
der to establish a relation between SSB of the rotating 
mean field and an appearance of the corresponding RPA 
spurious mode (cf [13). While a qualitative discussion 
about a relation between SSB effects and RPA spurious 
modes in nonrotating nuclei could be found in literature 
(cf ReffllH, here we present the first realistic quantita- 
tive attempt to get a thorough insight into this relation 
in rotating nuclei. The self-consistent determination of 
the constants G and K2 will be discussed below. 

Using the generalized Bogoliubov transformation for 
quasiparticles (for example, for the positive signature 
quasiparticle we have of = J^k^kiC^ + ^u'^k) ^-nd the 
variational principle (see details in Ref0) , we obtain the 
Hartree-Bogoliubov (HB) equations. The HB equations 
are solved using the " experimental" values of deformation 
parameters /? and 7 (see Fignj, obtained in Refill as an 
input. These values are extracted from experimental data 
[13. For the pairing field = -SkiGr{Pr) = 6kiAr{^) 
we assume a phenomenological dependence of the pairing 
gap At upon the rotational frequency Q 



Ar(o)[i-i(f )2] foTn<n, 

A,(0)i(^)2 ' for n>Q, 



(5) 



introduced in ReflU Here, Qc = 0.32, 0.33 MeV is a ro- 
tational frequency where the first band crossing (which 
is approximately the same for protons and neutrons) oc- 
curs for ^'"^^Dy and -'^^^Er, respectively. The values of 



pairing gaps A„(0) = 0.857 MeV, Ap(0) = 0.879 MeV 
(for i56Dy) and A„(0) = 0.874 MeV, Ap(0) = 0.884 MeV 
(for ^^^Er) are obtained from the odd-even mass differ- 
ence (see also Ref ITsi) ) . The Nilsson parameters n and /i 
are taken from systematic analyses for all deformed nu- 
clei .8i|. In our calculations we include all shells up to 
= 8 and consider AA^ = 2 mixing. The configura- 
tion space exhausted 97% of energy weighted sum rule 
for quadrupole transitions. 
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FIG. 1: (Color online) Equilibrium deformations in /3-7 plane 
as a function of the angular momentum / — (Jx) (in units of 
h). 

It is enough to solve the HB equations only for the 
positive signature quasiparticle energies since the neg- 
ative signature eigenvalues and eigenvectors (Wj,^, V^^) 
are obtained from the positive ones. Thus, our quasi- 
particle operators are defined with respect to the yrast 
state I) (the lowest HB state at given fi) which is our 
quasiparticle vacuum ai\) = 0, a-^l) = 0. 

In the limit of the harmonic oscillator, at the self- 
consistent energy minimum the expectation values of the 
operators Qo and Q^"^ expressed in the double-stretched 
coordinates, Eqs.l^j), are zeros 



(Qo) = (Q^+^)=o 



(6) 



In other words, the self-consistent residual interaction 
does not change the mean field equilibrium deformation. 
To check the self-consistency of the HB solutions, we cal- 
culate Eqs.® when all terms are included in the HB 
equations. First, we would like to stress that the "experi- 
mental" values of /3 and 7 correspond, indeed, to the min- 
imum of the total mean field (HB) energy. The variation 
of the deformation parameters (3 and 7 around the "ex- 
perimental" (equilibrium) values results in the increase 
of the HB energy. This minimum becomes very shallow, 
however, with the increase of the rotational frequency. 
Second, double-stretched quadrupole moments are ap- 
proximately zero for all values of the equilibrium defor- 
mation parameters. A small deviation from the equi- 
librium deformation leads to a strong deviation of the 
double-stretched quadrupole moments from zero values 
and correlates with the increase of the HB energy. The 
main outcome of the HB calculations is the onset of the 
non-axial deformation around Tifi ~ 0.32 MeV where the 
backbending takes place in both nucleus. While our re- 
sults are similar to the prediction of Ref 9 for A^ ^ 90, we 
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reproduce better the experimental data. Our equilibrium 
deformations are a result of the self-consistent solutions 
of the HB equations, whereas the authors of Ref0 used 
fixed phenomenological inertial parameters. 

To describe quantum oscillations around mean field 
solutions the boson-like operators b'^^- = ct^'^t^^ki = 
a'^a^jb'^j = cej^af are used. The first equality in- 
troduces the positive signature boson, while the other 
two determine the negative signature ones. These two- 
quasiparticle operators are treated in the quasi-boson ap- 
proximation (QBA) as an elementary bosons, i.e., all 
commutators between them are approximated by their 
expectation values with the uncorrelated HB vacuum Q . 
The corresponding commutation relations can be found 
in Ref .& In this approximation any single-particle op- 
erator F can be expressed as = (F) + F^^^ + F^^^ 
where the second and third terms are linear and bilinear 
order terms in the boson expansion. We recall that in 
the QBA one includes all second order terms into the 
boson Hamiltonian such that {F - (F))^ = ^(1)^(1). 
The positive and negative signature boson spaces are not 
mixed, since the corresponding operators commute and 
Hn=Hn{a = +) + Hn{cT^-). 

The CRPA Hamiltonian is diagonalized by solving the 
equations of motion for each signature separately. As a 
result, we obtain the following determinant of the secular 
equations 

.F(c.A)=det (R-1 ), (7) 

which is the fifth and the second order for the positive and 
negative signature, respectively, and c = K2 or Gt- The 
matrix elements Rkm{uJ\) = Y.^, <lk.f^qm,t,Cfj^/{Ef^-ujl) 
involve the coefficients C^"' — E^ or lox for differ- 
ent combinations of matrix elements qk,ii (see details in 
Refsl6l ll7i) . The zeros of the function J-{lox) = yield the 
CRPA eigenfrequencies lox. Since the mean field violates 
the rotational invariance and particle number conserva- 
tion law, among the CRPA eigenfrequencies there exist 
few spurious solutions. In this paper we focus our atten- 
tion on the SSB effects related to the rotation, since the 
pairing can be treated in the same way. 

Introducing the operator r+ — (Jy^^ — iji^^)/ \/2{J^ 
such that 

[Hnia^-),T+] = nT+, [r,r+] = i, r = (r+)+ 

. 

one is able to separate negative signature vibrational 
modes from the "spurious" solution at ujx = f2. Equa- 
tions © describe a collective phonon that creates a col- 
lective rotation. This phonon is related to the symmetry 
broken by the external rotational field (the cranking term 
in Eq.^). 

The other spurious solutions are associated with the 
rotation around the x axis and the particle number con- 
servation law, 

[Hn,X] = [Hn,Nr] = 0. (9) 



The mode associated with the rotation about the x axis 
allows one to determine the Thouless-Valatin moment of 
inertia JTtv using the positive signature term of the full 
Hamiltonian 

[Hn{a^+),i^^^, (10) 

JTV 

Here the angle operator $ is the canonical partner of the 
angular momentum operator Jx- A similar procedure 
can be applied for the second spurious mode in Eqs.0 
to obtain the mass parameters for neutron or protons jl8L 

It is important to hold a self-consistency at the CRPA 
level as well as at the mean field. In the harmonic pscilla- 
tor limit the self-consistent constant, K2 = ^-^j, war- 
rants the fulfillment all conservation laws in the CRPA 
for rotating nuclei Therefore, in our realistic calcula- 
tions we define the constants from the requirement of the 
fulfillment of the conservation laws, Eqs.®, and a sepa- 
ration of the rotational mode from the vibrational ones, 
Eqs.llHl). It was proved in Refill for the self-consistent 
model, which can be solved exactly, that the dynamical 

(2) 

Jfjg and the Thouless-Valatin Jtv moments of inertia 
must coincide, if one found a self-consistent mean field 
minimum and spurious solutions are separated from the 
physical ones. Our results (see Figl^J demonstrate a good 
self-consistency between the mean field and the CRPA 
calculations, indeed. 
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FIG. 2: The rotational dependence of the dynamical J^jg — 
—d^EHs/dQ^ — d{Jx)/dft (a dash line) and the Thouless- 
Valatin Jtv moments of inertia (a solid line). Here, Ehb is 
a mean field value of the full Hamiltonian. 

To analyze the low-lying excited states we construct 
the Routhian function for each rotational band {ly — 
yrast, /3, 7, .... ) R^{n) = E^{n) - ni^{n) and de- 
fine the experimental excitation energy in the rotating 
frame u;,y{il.) = Ru{^) — Ryr{^) [l^l- This energy can 
be directly compared with the corresponding solutions 
hujx of the CRPA secular equations. The experimental 
Routhians Ri, are obtained from experimental data (211] . 
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The results demonstrate a good agreement between the- 
ory and experiment (see Fig. 3). The lowest collective 
7 vibrational frequency for the positive signature states 
(even spins) becomes zero at hftcr ^ 0.324 and 0.33 MeV 
for ^^^Dy and ^^'^Er, respectively. As discussed above, 



persists. 



near the rotational frequency flc 
curs in the considered cases. 



the backbending oc- 
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FIG. 3: (Color online) The excitation energies in the rotating 
frame tjy(f2) for the positive signature 7 (top panel) and P 
vibrational (bottom panel) states. The results of calculations 
are connected by a solid line. The lowest two-quasiparticle 
states shown at the upper panel are connected by a dash line. 
The experimental values for different excitations are indicated 
by open circles connected by a thin line to guide eyes. 

In order to understand this correlation, let us consider 
an axially deformed system, defined by the Hamiltonian 
H in the laboratory frame, that rotates about a symme- 
try axis z with a rotational frequency fl. The angular mo- 
mentum is a good quantum number and, consequently, 
[Jz,0]^] = K0\^. Here, the CRPA phonon O]^ describes 
the vibrational state with K being the value of the an- 
gular momentum carried by the phonons along the 
symmetry axis, z axis. Thus, one obtains 



[Hn, 0\.] = [H- ah, 0|^] = {CuK - KVt)O^K = ^kO^- 

(11) 

This equation implies that at the rotational frequency 
ilcr — i^k/ K one of the CRPA frequency ljk vanishes in 
the rotating frame (see discussion in Refs l2^l23j) . At this 
frequency we could expect the SSB effect of the rotating 
mean field due to the appearance of the Goldstone boson 
related to the multipole-multipolc forces with quantum 
number K. For an axially deformed system one obtains 
the breaking of the axial symmetry, since the lowest crit- 
ical frequency corresponds to 7 vibrations with K = 2 
I22I l2a| . In contrast to nonrotating case [3, the RPA 
does not break down at the bifurcation point, since in 
the laboratory system the corresponding vibration still 
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FIG. 4: The rotational dependence of the mean field value 
(Jz) (in units h) for ^'''^Dy (left panel) and ^^**Er (right panel). 
The noncollective angular momentum is built by the align- 
ment of the separate two-quasiparticle pairs till Qci ~ 0.32 
(MeV) in both nuclei. 

Guided by this analysis, we solve the HB equations 
with rotation about the symmetry axis z (the noncollec- 
tive regime). Since the angular momentum of a quasi- 
particle state is conserved, the lowest two-quasiparticle 
configuration with the largest deformation aligned or- 
bitals builds the total angular momentum of the sys- 
tem when it crosses the ground state configuration with 
zero angular momentum. The total angular momen- 
tum increases each time stepwise by an amount of an- 
gular momentum carried by a two-quasiparticle config- 
uration that crosses the renewed ground state configu- 
ration. For ^^^Dy (see FigQ} (Jz) is zero for all rota- 
tional frequencies in the range of values < Tifl < 0.275 
MeV. At nn = 0.275 MeV the energy of the neu- 
tron two-quasiparticle state originating from the shells 
^9/2 ^ /7/2 (for = A = it is the configuration 
built from the Nilsson states |[505] (g) |[514]) goes to 
zero. The nuclear angular momentum is determined by 
the angular momentum of this two-quasiparticle state: 
{Jz) = I + i = Note that one would expect the 
onset of the nonaxiallity at this frequency, according to 
the quasiparticle picture. According to our results, at 
hil = 0.31 — 0.32 MeV the expectation value (Jz) is in- 
creased stepwise by next four units. The system remains 
an axially deformed creating a state with (Jz) = I2h. 
At hilcr = 0.324 MeV the system acquires additional 
two units and the shape transition occurs to a nonaxially 
deformed shape. We recall that at this frequency 7 vibra- 
tions of the positive signature carrying two units of the 
angular momentum vanish. Comparing the calculated 
value of the critical frequency with the experimental es- 
timation hn„ = TiCjk=2/2 = 0.691/2MeI^ - 0.345 MeV, 
one could find a good agreement between theory and ex- 
periment. As was mentioned above, the potential energy 
surface is very shallow at large rotational frequencies. In 
fact, the energy minima for a rotation around the axis 
z and the axis x are almost degenerate. The difference 
is about 15 keV near the critical rotational frequency. 
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At the bifurcation point a competition between the col- 
lective (around the axis x) and noncoUective rotations 
breaks the axial symmetry and leads to nonaxial shapes. 

For ^^^Er the experimental critical value is expected 
at hn„ = hCjK=2/K = 0.613/2 MeV - 0.306 MeV. 
The angular momentum is zero up to the rotational fre- 
quency Tifl — 0.272 MeV. At this frequency the proton 
two-quasiparticle state originating from 177 /2 ® /2 shells 
(for = A = this configuration is built from the Nils- 
son states I [404] (g) |[402]) contributes to the value of the 
angular momentum | -I- | = 6h. The system remains the 
axially deformed, getting another two units of the angu- 
lar momentum, up to hflcr ~ 0.33 MeV, where 7 vibra- 
tions of the positive signature are vanished in the rotating 
frame. Again, there is a tiny energy difference between 
the collective and noncoUective rotations (~ 15keV). At 
this frequency the expectation value (Jz) is increased by 
two more units and for Q > Qci the system is driven into 
the domain of triaxial shape. 

According to the CRPA 6], the reduced transition 
probability between 7 (one-phonon) states of the posi- 
tive signature B{E2,AI = 2) ~ |(Q^=+(A/ = 2))\^ is 
of the same order of magnitude as the collective electric 
quadrupole transitions between yrast (vacuum) states. 
In fact, in both nuclei the B{E2,AI = 2) strength ex- 
ceeds few tens of Weisskopf units. Detailed calculations 
of the transition probabilities are beyond the scope of the 
present paper and will be presented elsewhere. 

Summarizing, for the first time the cranking HB and 
RPA equations are solved in a self-consistent manner. 
We obtain a good agreement with available experimental 



data for low-lying vibrational excitations at high spins 
in ^^^Dy and ^^^Er. According to our analysis, the 
alignment decreases the pairing correlations and, con- 
sequently, the 7 vibrations are softening in axially de- 
formed nuclei. We stress that the two-quasiparticle align- 
ment does not create the backbending. Rotating around 
the axis which is perpendicular to the symmetry axis, 
we found that in ^^^Dy and ^^^Er the backbending oc- 
curs at the critical rotational frequency Qcr ~ ^k=iI'2 
where Cjk=2 is a 7 vibrational excitations in the labo- 
ratory frame at 17 = 0. At this frequency the positive 
signature 7 vibrational excitations vanish in the rotating 
frame, i.e., ljJk=2 = 0. As a result, the nuclear mean 
field spontaneously breaks the axial symmetry and gives 
rise to nonaxially deformed shape in the rotating frame. 
Thus, the inclusion of quantum fluctuations around the 
mean field approach extending our understanding of the 
backbending phenomenon. 
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